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Abstract. We present a simple model to account for the rheological behavior 
observed in recent experiments on micellar gels. The model combines attachment- 
detachment kinetics with stretching due to shear, and shows well-defined jammed and 
flowing states. The large deviation function (LDF) for the coarse-grained velocity 
becomes increasingly non-quadratic (and the probability distribution function non- 
Gaussian) as the applied force F is increased, in a range near the yield threshold. The 
power fluctuations are found to obey a steady-state fluctuation relation (FR) at small 
F. However, the FR is violated when F is near the transition from the flowing to 
the jammed state although the LDF still exists; the antisymmetric part of the LDF is 
found to be nonlinear in its argument. Our approach suggests that large fluctuations 
and motion in a direction opposite to an imposed force are likely to occur in a wider 
class of systems near yielding. 
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Fluctuation relations (FR) as originally formulated [H El El IH El El E] are exact 
statements connecting the relative probabilities of observing the production and the 
consumption of entropy at a given rate in a driven thermal system, and can be expressed 
as symmetry properties of the large-deviation function of the entropy production rate. 
An observable X^-, for example, the average power delivered over a time interval r is said 
to have the large-deviation property if its probability density decreases exponentially 
for large r, Prob(X^^oo = a) ~ exp[— rVV(a)]. The decay rate W is called the large- 
deviation function (LDF), and its behaviour for large a encodes information about 
the statistics of extremes of the underlying random process [H [9]. The steady-state 
fluctuation relation [21 13] states that >V(a) — >V(— a) oc a. 

In this paper we study numerically in some detail the large-deviation behaviour 
of a model for a macroscopic degree of freedom driven through a medium of dynamic 
attachment points. The model is motivated by precision creep rheometry studies [lO] 
of a micellar gel at controlled stresses below its nominal yield point, that revealed 
that the small positive mean shear-rate, i.e., in the direction favoured by the imposed 
stress, was composed of a highly irregular time-series of positive and negative shearing 
events. Moreover the shear-rate fluctuations, which at constant stress are also power 
fluctuations, obeyed [10] a fluctuation relation of the Gallavotti-Cohen type [21 [3] . The 
measurements reported were made on a macroscopic degree of freedom, the angular 
position of the rheometer plate, which could not be influenced perceptibly by thermal 
noise. Fluctuations in this system must be a consequence of the imposed drive. The 
probability distribution of the angular-velocity fluctuations was found to be strongly 
non-Gaussian at large imposed stress. An effective temperature extracted from a 
comparison to the Gallavotti-Cohen relation was found to increase with the applied 
stress. 

Our theoretical model reproduces the findings of the experiment but also finds 
departures from the fluctuation relation in the strict sense, in a certain parameter range, 
despite the existence of LDF. That is, it finds that the antisymmetric part of the LDF 
departs from linearity in its argument. Moreover, the model we present should apply to 
a wider class physical problems involving yield or escape in the presence of fluctuations. 
We shall return to these points at the end of the paper. 

Consider an external force acting on some coordinate of a thermal system with 
inverse temperature /3, and let the random variable Wt be the instantaneous rate of 
doing work in a given realization of the dynamics. Define 



r JO 

to be the rate of doing work, binned or averaged over a time scale r. The steady-state 
fluctuation relation (FR) [TH [T2] tells us 



where P{Wr) is the probability distribution function (PDF) of Wr, and x denotes 
asymptotic equality. Underlying (|2|) is a more general relation, the existence of the 
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large deviation function (LDF) [5]. The principle of large deviation consists in the 
existence of the limit 

lim --lnP(l^,) = W, (3) 

r— >-oo f 

where the limiting quantity W is called the rate function or large-deviation function 
(LDF). Eq. [2] is then a statement about a symmetry property of the LDF, viz., that its 
antisymmetric part is linear in its argument. A large class of experimental [T3| [El [15] 
and simulation [121 [IT] results and theoretical calculations [T^ f2U\ [2T] confirm the 
existence of relations of the FR type even in systems driven far away from equilibrium. 
Note that some of these systems are not obviously characterized by a thermodynamic 
temperature. In such cases the existence of a relation like (|2|) offers one way of defining 
an effective temperature [22] . 

The exact linearity in Wr as required by ([2]) is a strong restriction for a non- 
Gaussian random variable. Naturally there are situations [231 1211 [251 126] where it does 
not hold. In the example of [231 [2l] , however, the large-deviation property itself does 
not hold. In our simple model, by contrast, the conventional FR fails to hold in some 
regions of parameter space even though the LDF exists. The model has the additional 
virtue of being constructed to model a physical situation, rather than as a mathematical 
counterexample, and will therefore be of wide interest. 

Here are the main results of this work: (1) the model shows a sharp crossover from 
a creeping jammed state to steady flow (Fig. [2]); (2) deep in the jammed state, the 
velocity fluctuations are Gaussian and obey the FR (Fig. [3|); (3) near the threshold to 
free flow, the velocity fluctuations become non-Gaussian, but still obey the FR (Fig. [5]); 
all these findings are in conformity with experiments [10] . (4) Just below the threshold, 
the velocity fluctuations are non-Gaussian and they don't obey FR, even though the 
LDF exists as is shown in Fig. [1]). (5) If we keep the number of attachment points fixed, 
the effective temperature {Te/f) decreases with increasing applied force (F) in contrast 
to what is found in the experiment. To produce the correct T^ff vs F trend, we must 
allow the number of attachment points to increase with increasing F (Fig. [7]). 

The paper is organised as follows: in Sec [D we present our model and we discuss 
the mean-field version of the model in Sec. [2l We present details of our exploration of 
parameter space and the results in Sec. [3] and conclude the paper with a discussion in 
Sec. m 

1. The model 

The micellar gel sample in the experiment is taken in a rheometer with a cone-plate 
geometry to ensure uniform strain rate through-out the sample, and the upper plate 
is rotated with a constant torque while keeping the lower plate fixed. The detailed 
geometry of the rheometer, however, is not important for the observed findings. We 
model the micellar gel medium as a collection of springs that are stretched by the 
applied torque when they are attached to the plates, and can detach when stretched by 
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a high enough force. We do not associate the "springs" with individual molecules or 
micelles but rather with adhering, deformable domains in the material, whose size we 
do not know. We assume the springs always remain attached to the stationary plate. 



1 If^ r^^r^i^ 



Effective 1-D geometr 




Attached 
Detached 



Figure 1. A schematic illustration of the model. The "springs" (sec text) can get 
attached to or detached from the upper plate. The experimental geometry can be 
thought of as effectively one-dimcnsional with the force being applied in a particular 
direction. 



We use an effective one-dimensional description in which X{t) is the total (angular) 
displacement of the upper plate as in Fig. [Hand let V{t) = dX{t) / dt be its instantaneous 
velocity. Let us consider fi{t) = kiXi be the force on the zth spring at time t where ki 
and Xi are respectively the spring constant and the extension of the spring. Then the 
spring will pull the plate backwards only if it is attached with the plate. Thus, we can 
write down the equation of motion for the upper plate as 

.dV{t) 



M- 



dt 



+ BVit) = F-J2s,m, 



(4) 



where M is the mass of the plate, B a viscous damping coefficient, F the force (actually 
torque) on the rheometer plates and Sj is a two state variable which can take on values 
or 1. If Sj = 1, then the ith spring is attached with the plate and it is detached 
otherwise. Note that we have not specified the number N of "springs" in the model. 
It is not clear how to do this in the absence of a detailed microscopic theory. It is 
entirely possible that depends on the imposed torque, or even that it is determined 
dynamically. We will assume it is a parameter of the system, and show the behavior in 
the N — F plane. 

Now, if the spring is attached to the plate, it will stretch with the velocity of the 
plate; if it is detached it will relax. Thus 

dxiit) , 



dt 



Si)^kiXi{t) + SiVit) 



where 7 is a kinetic coefficient. Thus, 

dm) 



dt 



{l-s,)^kJ,{t) + SikiV{t). 



(5) 



(6) 
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Assume for simplicity that all /cj's are equal, define 7/cj = 1/r, redefine kiV V in Eq. 
(|6]) and B/ki ^ B in Eq. (jl]). Ignoring inertia, we then obtain equations for V and 

BV{t) = F-Y.s,m) (7) 

i 

^ = -{l-s,)m/r + s,V{t). (8) 

One could imagine more complicated modes of relaxation, for example the springs can 
partly redistribute forces among themselves. But Eq. ^ is the simplest possible model 
that contains the dominant mechanism. We will discuss the effect of a diffusive term 
later. 

The state variables Sj are assumed to follow a stochastic dynamics. Let 
Pi{t) =Prob(sj = be the probability that the ith spring is attached at time t. 
Then 

= -WDP^{t) + Wa{1 - m)), (9) 

where Wa and Wd are the attachment and detachment rates, the most important input 
parameters of the model. It is possible to engineer these rates to reproduce different 
behaviours by the model. The feature that is essential to get negative fluctuations is 
that a spring with a large force on it gets reattached and pulls the plate in the opposite 
direction. 

The fluid is a highly viscous glassy material and we can safely ignore the inertia 
term in Eq. (jlj). For notational convenience we set B to unity so that V = F — Sifi- 
For a suitable choice of parameter values, as will be shown below, the model shows a 
jammed- flowing transition. This is not a true phase transition, but a strong crossover 
from slow creep to free flow. 

2. The mean field calculation 

In the mean field approximation, we assume the forces on all springs are the same and 
that they take some mean value in the steady state, as do the state variables (sj). In 
this approximation, (/j) = / and (sj) = s so that 

{V) = F- Y.{f,s,) ^ F - Y.U^){s.) = F- Nsf, (10) 

i i 

where is the total number of springs. In the steady state, ([8]) will yield 

(1 - s)f/r = s{V) ^ / = ns{V)/il - s). (11) 
Using the above relation in Eq. ( |TOl) . we find the velocity 



(K) . -ii-^. (12) 

1 — s) + As^r 



Yielding and large deviations in micellar gels: a model 



6 



In two extreme limits, if s = 1, {V) = and if s = 0, (V^) = F. In the steady state, 
dPi{t)/dt = so that 

Therefore, the steady state value of s, within the mean-field approximation, will be 

s = Is,) = y s,P, = Pi = — . (14) 

The mean force on each spring becomes 

f = W , ■ (15) 

(1 — s) + Ns^T 

Thus, we see that {V) depends on three time scales. To simplify the discussion, we 
fix r and Wa and take Wd/Wa to have an activated form. We recall that the mechanism 
of having a negative velocity (in the direction opposite to F) events is that a spring with 
a large force on it gets reattached to the plate before it has completely relaxed its force. 
This can happen if the force relaxation is much slower than the attachment-detachment 
kinetics of the springs. Thus, to have a large number of negative events, we must have 
tWa ^ 1. To ensure this, we chose the parameters as follows: r = 2.5, Wa = 100 and 
Wf) = Wag"'^-^^' with a = 2.0 and /o = 1.0. The springs can get attached with the 
plate at a constant rate irrespective of the force on it. But if it is already attached to 
the plate, it is more probable to get detached as the force on the spring increases. We 
have introduced the parameter a to obtain a reasonably sharp transition from jammed 
to flowing state, at a force whose value is controlled by /q. Restrictig the form of Wd 
allows us to worry about one less parameter of the model. The predictions of mean-field 
version of the model are shown with solid lines in Fig. |2l 

The activated nature of the attachment-detachment is reminiscent of rheology 
models with traps such as the SGR [271 EHl EH]. An essential difference between our 
model and the trap models [271 128] is that the springs in our model retain the forces on 
them even after getting detached from the plate whereas the strain on a spring in the 
trap model becomes zero after it comes out of a trap. The simplest trap models can not 
show negative velocity events. 



3. Details of the simulation and the results 



We simulate the equations ([7])- ([9]), setting B = 1, through the kinetic Monte-Carlo 
(KMC) method [SOj to obtain the behaviour of our model. The advantage of KMC over 
the conventional Monte-Carlo method is that the time scale of the dynamics is entirely 
determined by the various rates of the problem. As we have discussed in the previous 
section, the various parameters of the model were chosen using the mean-field analysis 
to locate a regime in which a clear transition from jammed to flowing was seen, and 
adjusting rates to allow for significant negative velocity fluctuations. 

The model shows a well-defined jammed to flowing transition that becomes sharper 
as we increase the number of springs N as is shown in Fig. |2] where the symbols are the 
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Figure 2. (a) The flow curves obtained from tlie simulation for various number 
of springs with a particular set of parameter values (see text), the lines are the 
corresponding mean-field solution, (b) The flow curve is shown for TV = 200 for clarity. 
The mean-field solution overestimates the threshold. The regions are roughly marked 
based on whether the fluctuations are Gaussian (G) or non-Gaussian (NG) and whether 
they do obey the conventional FR or not. 



simulation values and the corresponding curves are the mean-field solution. In Fig. [2](b) 
we show the fiow curve for = 200 and the inset shows the behaviour in the jammed 
state in a semi-log plot where it is evident that even in the jammed state, the velocity is 
actually non-zero. The mean-field solution overestimates the threshold but agrees well 
with simulation results away from the transition. 

In the jammed state, there is a significant number of negative velocity events, 
however, there are none once the system goes to the fiowing state. To understand the 
fiuctuations in the jammed state and test the regime of validity of the FR, we keep 
= 200 fixed, and find a threshold force F = 97.0. We take a value of the external 
force F = 20.0 deep in the jammed state. From the instantaneous velocity as shown in 
Fig. [3] (a), we see that there is a significant number of negative velocity events. Since the 
applied force is constant, the statistics of velocity and power fiuctuations are the same. 
We denote the velocity (or power) fiuctuations with respect to the mean, averaged over 
a time interval r by Wr'- 

1 /■'^+^ V(t) 

(16) 

where the initial time T is arbitrary but taken in such a way that the difference between 
two consecutive T is greater than the mean correlation time of velocity. The PDF 
of Wr is shown in Fig. [3]^b). The overlap of curves in the scaled plot [Fig. |3]^c)] is 
evidence of the existence of the large-deviation function for the velocity. The nature of 
the distribution is determined by two parameters, Skewness (S) and Kurtosis {K): 

S = {6Wr)ycT^ and K = {SWry/a^ - 3, (17) 

where SWr = Wr — (Wr) and cr^ = {6W^). S and K are defined so as to be zero for a 
Gaussian process and Fig. [3]^d) shows that the distribution is Gaussian. Also, Wr obeys 




Figure 3. The behaviour of the model for the particular set of parameter values as 
specified in the text with N = 200 and F — 20.0 (a) The instantaneous velocity as a 
function of time, the thick line denotes the average velocity. There are a significant 
number of negative velocity events, (b) The probability distribution function (PDF) 
of the scaled coarse-grained velocity or work fluctuation (both are same since the 
applied force is constant) Wr- (c) The large-deviation function (LDF) for Wr- (d) The 
skewness (S) and Kurtosis {K) for velocity fluctuation are very small implying the PDF 
is Gaussian, (e) \n[P{Wr)/ P{-Wr)] vs Wr for various r. (f) (1/r) ln[P(W^)/P(-W^)] 
vs Wt for various r collapse to a master curve that is straight line signifies that the 
velocity fiuctuation obeys fluctuation relation. 



the fluctuation relation as is evident from the plot of ln[P{Wr)/P{—Wr)]/T vs Wr for 
various r collapsing to a master curve that is a straight line going through the origin 
[Fig. Hf)]. 

One of the interesting features of the model is that if we are very close to the 
threshold, even though the large deviation function exists, Wr doesn't obey the standard 
fluctuation relation, as is seen in Fig. HI Here we keep the applied force F = 96.0 which 
is very close to the threshold value. In this case, the velocity fluctuations become non- 
Gaussian as is evident from the large values of Skewness (5*) and Kurtosis (K) in Fig. 
m^b). The plot of \n[P(Wr) / P{—Wr)] vs Wr deviates from a straight line and if we 
scale ln[P{Wr) / P{—Wr)] by r, even though we obtain data collapse, the master curve 
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Figure 4. The number of springs is iV = 200 and F = 96.0 which is very close to 
the threshold value, (a) The large-deviation function exists for the velocity fluctuation 
Wr- (b) The large value of S and K imply the PDF for WV is non-Gaussian, (c) 
hi[P{Wr) / P{—Wr)] vs Wr for various r as shown in the figure. The curves deviate 
from the expected straight line if FR was obeyed by Wr- (d) When we scale the various 
curves in (c) by r, they show excellent data collapse, but the collapsed data deviates 
significantly from a straight line implying the violation of the conventional FR. The 
straight line in the figure is just a guide to the eye. 

is no longer a straight line as shown in Fig. Wid). As long as we are very close to the 
threshold force, similar nonlinear FR curves are obtained even if we change the number 
of springs A^. Our simple model thus offers an example of a system with substantial 
negative fluctuations and excellent data collapse consistent with the large-deviation 
property, but in the which the antisymmetric part of the large-deviation function is 
strongly nonlinear. 

However, if we move slightly away from the threshold but still within the jammed 
state, the velocity fluctuation remains non-Gaussian, but the fluctuation relation is 
obeyed over the entire range of our data (Fig. [5]). This shows that the source of 
deviation from the fluctuation relations is not merely the non-Gaussian nature of the 
velocity fluctuations, but it is an intrinsic nature of the model and stems from a complex 
mechanism near the threshold. We do not completely understand the origin of such 
behaviour, but further work in this direction should elucidate this very interesting 
phenomenon. 




Figure 5. The number of springs is = 200 and applied force F = 86 which is close 
to the threshold value, (a) The large deviation function exists, (b) The non-zero value 
of S and K imply that the PDF of velocity fluctuation is non-Gaussian, (c) Plot of 
ln[P{Wr) / P{—Wr)] vs Wr for various t as shown in the figure, (d) When we scale 
ln[P{Wr) / P{—Wr)] with T, the curves show data collapse and the master curve is a 
straight line going through the origin ascertaining the validity of FR. 

NG+FR NG+NFR 

U No negative 

Deep in the jammed state Threshold ^ 



Figure 6. For a certain number of springs, the system goes from a jammed creeping 
state at low external force (F) to a free flowing regime as F is increased. Deep in 
the jammed state, the PDF of velocity fluctuation is Gaussian (G) and they obey 
FR. However, as we increase F towards the transition, close to the transition, PDF 
becomes non-Gaussian (NG) although FR is still obeyed. Very close to the threshold, 
the PDF is NG and FR doesn't hold anymore. We do not see any negative velocity 
events beyond the threshold. 



Let us summarize the behaviour of the model in a schematic phase diagram, Fig. 
ini Deep in the jammed state, the velocity fluctuations are Gaussian (G) and they 
obey the fluctuation relations (FR). Near the threshold force of the jammed-to-flowing 
transition, the velocity fluctuations become non-Gaussian (NG) but they still obey FR. 
As we approach the transition, very close to the threshold, the non-Gaussian nature of 
the velocity fluctuation remains but they no longer obey FR. We note that these features 
of our model are similar to the statistical properties of entropy-consuming fluctuations 



Yielding and large deviations in micellar gels: a model 



11 




F 



100 



Figure 7. (a) The effective temperature Teff extracted from Eq. ([2]) in arbitrary 
unit is plotted as a function of N and F. If we keep N fixed in the simulation, T^^ff 
decreases with increasing F as is found from the two paths shown in the figure by two 
blue solid arrows corresponding to = 200 and N = 400. However, if we consider 
a different path, one such possible path is shown by the red dotted arrow, where A'' 
increases with increasing F, Teff increases on this path, (b) Teff vs F for the path 
denoted by the red dotted arrow in (a) is shown for clarity. The number of filaments 
corresponding to a particular F is listed in Table [TJ 



in jammed states of laponite suspensions [25] . 



The slope of the scaled FR plot, analogous to (3 in Eq. (|2]), can be thought of as 
the inverse of an effective temperature. In the experiment, the effective temperature 
[Teff) increases as F increases [10]. We have pointed out earlier that there is no reason 
for the number of attachment points N to remain fixed in the model. In fact, it is more 
reasonable that changes with F, since large applied force can disentangle or break 
micelles or disrupt domains giving rise to more independent regions in the dynamics. To 
see the behaviour of T^ff as a function of F within our model, we have plotted T^ff as a 
function of both F and (Fig. [7]). Let us first see what happens if we keep N fixed. As 
we have shown in Fig. [Tj^a) for A^ = 200 and A^ = 400 (two solid arrows), T^ff decreases 
as F increases with constant A^. This is in complete contrast to what was found in the 
experiment [10]. However, as A^ increases, T^ff increases. Thus, to be consistent with 
our model, it must be that the system moves on a path on which A^ changes with F. 
With this in mind, it is possible to identify a path in the 3d space of (A^, F, T^ff) where 
the model reproduces the correct trend of Teff vs F. One such possible path is shown 
by the dotted arrows in Fig. [Tj^a) and the corresponding T^ff vs F behaviours is shown 
in Fig. O^b) for clarity. We also show in Table [1] the particular number of attachment 
points for a particular force F corresponding to this path. We emphasize that this is 
not the only possible path consistent with the experimental trends, and the particular 
path that the experiment will follow is going to depend on the microscopic details of 
the experiment. 
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F 20 40 60 80 100 
N 50 100 150 400 700 



Table 1. The number of springs N taken at a particular F 



4. Discussion and conclusion 

In this work we have presented a simple model to understand a particular set of 
experiments where it was found that the velocity (or power, since the applied force is 
constant) fluctuations obey the fluctuation relations. The force- dependent attachment- 
detachment kinetics of the springs with the plate is the main mechanism behind the 
observed negative velocity. When the applied force is very close to the threshold, the 
probability distribution function of velocity fluctuations becomes non-Gaussian and the 
strong departures from a conventional fluctuation relation are seen. This is especially 
interesting given that the large-deviation property continues to hold, and we obtain data 
collapse when we plot \rL[P(Wr)/ P{—Wr)]/T as a function of Wr, though the master 
curve is not a straight line. To our knowledge, our model is the first physical example 
of such a violation of FRs where the LDF exists. We need more experiments and 
theoretical analysis to understand the origin of such a phenomenon. The observation of 
a linear dependence of\og[P(Wr)/P{—Wr)] over a range of Wr is not in itself our major 
finding, as a function that goes through zero will normally have a linear range. It is of 
greater significance (a) that the function goes through zero and has appreciable weight 
at negative arguments, suggesting that the model captures some of the essential physics 
of the experiment and (b) that we observe good data collapse even when the symmetry 
function departs from linearity. A worthwhile future direction will be to sample the rare 
events [SI |32] to improve statistics in the tail of the distribution, possibly elucidating 
the nature of this violation of FR. 

If we keep the number of springs taking part in the dynamics fixed, the observed 
trend of the effective temperature as a function of applied force is opposite to what was 
found in the experiment. However, it is more reasonable to vary the number of springs 
as F changes, since larger applied force may break entanglements, rupture micelles, or 
disrupt adhering domains. This allows the model to reproduce the correct trend oiT^ff. 

In the model, we have allowed a simple local relaxation mechanism for the springs. 
One can imagine more complicated modes of relaxation, for example, we can allow the 
springs to redistribute forces among their neighbors: 



We find that the presence of such a diffusive term doesn't affect the behaviour of the 
model much. During the attachment-detachment kinetics, other processes that could 
play a role include: spatial inhomogeneity and temporal variation of spring stiffness and 
their modification by local stretching and release, and the interplay of micellar lengths 
and relaxation time with imposed stresses [33] . 
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Finally, we expect that fluctuations near yielding in a wider class of systems could 
show features similar to those discussed here. We have in mind situations such as the 
dislocation- mediated flow of stressed crystals at non-zero temperature [3l] , the flow of 
glass through the mechanism proposed by Sausset et al. [35], and thermally assisted 
depinning in general [36] • Fluctuations that take a region from the downhill to the uphill 
side of a pinning barrier, which are clearly more likely to happen near yielding, where 
effective barriers are small, should give rise to negative-velocity events. 
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